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We construct model wavefunctions for the collective modes of fractional quantum Hall (FQH) 
systems. The wavefunctions are expressed in terms of symmetric polynomials characterized by a root 
partition that defines a "squeezed" basis, and show excellent agreement with exact diagonalization 
results for finite systems. In the long wavelength limit, we prove the model wavefunctions are 
identical to those predicted by the single-mode approximation, leading to intriguing interpretations 
of the collective modes from the perspective of the ground state guiding center metric. 



One of the main driving forces in characterizing the 
topological phases arising in the context of the fractional 
quantum Hall (FQH) effect jl[ has been the construction 
of the model wavefunctions and their parent Hamiltoni- 
ans ■ Depending on the nature of the effective inter- 
action within a partially filled Landau level, numerous 
phases appear, ranging from crystals to incompressible 
liquids with fractional Q, and possibly non-Abelian Q 
quasiparticle statistics. Within the lowest Landau level 
(LLL), the kinetic energy of the electrons is a constant, 
and all dynamics comes from the Coulomb repulsion . 
Thus the essential features of the many-body ground 
states and the low-lying excitations are solely determined 
by the particle statistics and quantum fluctuations that 
minimize the repulsion. Laughlin's wavefunction 0] is 
a prominent example: the statistics is taken care of by 
the odd power in the Jastrow factor, and for rotation- 
ally invariant system, the energy is minimized by placing 
all zeros of the wavefunction on the electrons. Subse- 
quently, an intuitive picture of "composite fermions" was 
put forward by Jain [6[ , that lead to a classification of the 
numerous fractional states observed in experiment, and 
provided a numerical procedure for constructing various 
model wavefunctions @. While much effort has been 
devoted to formulating FQH ground-state wavefunctions 
and those of charged excitations (quasihole and quasielec- 
tron), relatively little progress has been made in under- 
standing the neutral excitations since the seminal work 
of Girvin, MacDonald and Platzman that introduced the 
single mode approximation (SMA) to describe the low- 
est excitation in terms of a neutral density wave or the 
"magneto-roton" Very recently, this question was 
revisited and explicit wavefunctions for the neutral col- 
lective excitations were proposed for the Abelian, as well 
as non-Abelian, FQH states utilizing a multicomponcnt 
composite fermion approach 0, Q . 

In this Letter we construct model wavefunctions for the 
collective modes of FQH states by extending the Jack 
polynomial description [Hi EH]. We model the lowest 



neutral excitation as a dipole formed by a single quasi- 
electron and a single quasihole. At the Laughlin v = 1/3 



filling, we obtain a single bosonic mode that corresponds 
to the magneto-roton; in the Moore- Read v = 5/2 case, 
we obtain in addition a neutral fermionic mode fl3L ITij 
that stems from the non-Abelian statistics. We iden- 
tify the long-wavelength limit of the bosonic modes as 
a "spin-2" excitation (analogous to the graviton), and 
that of the fermionic mode as a "spin-3/2" (with a pos- 
sible analogy to the gravitino). Our model wavefunc- 
tions show excellent agreement with exact diagonaliza- 
tion results at all wavelengths. In the limit of long wave- 
lengths, we present a proof that the model wavefunctions 
become equivalent to the SMA result, which thus remains 
an accurate description even at energies lying above the 
threshold of the roton-pair continuum. 

Wc first review some basic properties of the Jack poly- 
nomials. The Jacks are symmetric multivariable polyno- 
mials Jx( z i, z 2, • • • z n c ), parametrized by a number a and 
a root partition A with length l\ < N e , N e being the num- 
ber of electrons. For bosonic FQH states, the statistics of 
the system restricts us to only symmetric monomials with 
a fixed total degree, which corresponds to a fixed angular 
momentum on a disk or a sphere [l5j . Fermionic FQH 
states are obtained by multiplying the bosonic counter- 
part with a Vandcrmonde determinant. The root par- 
titions encode the clustering properties that effectively 
"keep electrons apart" from each other; a "(fc,r, in- 
admissible" partition is such that no more than k elec- 
trons are found in r consecutive orbitals. The root par- 
tition tells us how the wavefunction vanishes when parti- 
cles are brought together. For example, the root partition 
of the ground state and the quasihole states at v = l/m 
Laughlin filling contains no more than one electron in 
m consecutive orbitals. The corresponding Jack vanishes 
at least with a power of m when two electrons approach 
each other. 

From a practical point of view, the Hilbert space we 
are interested in, after reduction by particle statistics and 
symmetry, only needs to include polynomials with appro- 
priate clustering properties (whereby Jacks are special 
cases). In fact, for the ground state and single quasihole 
states of the Read-Rezayi series (including the Laugh- 
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lin and Moore- Read states), such Hilbert space is one- 
dimensional and spanned by single Jack polynomial with 
a = — z r-j that can be obtained directly via recursive re- 
lations [L6|. Quasielectron states, on the other hand, are 
more complicated [l2[ because they contain local defects 
where electrons are forced to get closer to each other than 
allowed in the ground state. The same difficulty arises in 
collective modes which consist of quasielectron-quasihole 
pairs. In these cases the reduced Hibert space is no longer 
one-dimensional with a single Jack polynomial. 

We now proceed to construct explicit model wavcfunc- 
tions for the magneto-roton mode in the Laughlin state 
at v = 1/3 filling. We work in the spherical geom etry 
with a monopole of strength 2S placed at the center [15| . 
The root partition for the Laughlin ground state is well 
known: A = {100100100100 ••• 1001} [nj. The cluster- 
ing property in this case means that no more than a sin- 
gle electron can exist in three consecutive orbitals fl~7j ]. 
The collective mode in the limit of momentum k — > 
consists of two quasielectron-quasihole pairs, forming a 
quadrupole. As the momentum increases, a dipole mo- 
ment develops with the separation of one quasielectron- 
quasihole pair. This is summarized in the explicit set of 
root partitions as follows: 

1100001001001001001 • ■ -L = 2 
1100010001001001001 • • -L = 3 
1100010010001001001 • • -L = 4 
1100010010010001001 • • • L = 5 
1100010010010010001 • • • L = 6 

(1) 

We label the states by their total angular momentum 
L on the sphere. Note the ground state with the root 
partition {1001001 . . .1001} has L = 0, and the excita- 
tion with the smallest momentum that can be created is 
L = 2. In Eq.([T]) the black dot schematically indicates 
the position of a quasielectron, while the white dot that 
of a quasihole. To determine the position of a quasiparti- 
cle, we look at any three consecutive orbitals in the root 
partitions above, and count the number of electrons to 
see if it violates the ground state clustering property. In 
this particular case, if there is more (less) than one elec- 
tron, we then have a quasielectron (quasihole), which is 
located right below the middle of the three consecutive 
orbitals . Due to rotational invariance on the sphere, we 
next impose the highest weight condition on the wave- 
functions \4>\) to single out the state with quasiparticlcs 
piled up at the north pole: 

L + \tf) =0, 

where arc monomials with partition fi (Toi | . The 
summation is over all partitions \i that can be squeezed 



TABLE I. The overlap of the approximate model wavefunc- 
tions constructed from product rules and the true model wave- 
functions. 
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89.83% 


90.13% 


90.31% 


90.42% 


L=3 


86.42% 


86.99% 


87.37% 


87.63% 


L=4 


83.63% 


84.59% 


85.23% 


85.69% 



from the root partition A. For example, m{iooi} ~ 
z\ — z|i m {ono} ~ z i z 2 — Z\z\, and {0110} is squeezed 
from {1001}. The constraints in Eq. (|2j) substantially re- 
duce the Hilbert space dimension (e.g., the basis dimen- 
sion is less than 20 for 10 particles). The resulting lowest- 
energy eigenstates of the Hamiltonian, restricted to this 
Hilbert space, are very good approximations to the exact 
magneto-roton mode. 

The innovation we implement here is to impose an ad- 
ditional constraint that can be formally expressed as 

tW 2 |V>A> = 0. (3) 

Here V\ is the operator corresponding to the first Haldane 
pseudopotential 0] , and Cj annihilates an electron at zth 
orbital. This additional constraint renders \"4>x) unique 
by enforcing the following clustering property: the wave- 
function is vanishing only when two or more clusters of 
two particles coincides in real space. 

The resulting implementation is numerically much less 
expensive, with variational energies only slightly above 
the ones obtained in the Hilbert space defined by con- 
straints (2), and improving with the increase in system 
size. We note that the model wavefunctions \4>\) in- 
herit rich algebraic structures from the underlying Jacks. 
When the geometric normalization factors on the sphere 
are removed, the coefficients of the decomposition in Fock 
space are integers, with the coefficient of the root con- 
figuration normalized to 1. Furthermore, they satisfy a 
"product rule" 0, [3] if the first five orbitals are treated 
as one "big" orbital, which allows us to generate a large 
subset of coefficients recursively This suggests the prod- 
uct rule is not restricted to pure Jacks, and is robust 
against local defects of the wavefunction. An approxi- 
mation to can be built from the product rules; the 
overlap between the approximate and exact model wave- 
functions is high and increases with system size (see Ta- 
ble I). These properties reduce the computational cost 
for generating these model wavefunctions compared to 
direct diagonalization. 

We now proceed to evaluate the variational energies of 
the model wavefunctions obtained via constraints Eq.Q 
and j3]). In Fig. [TJ the variational energies are plotted 
versus momentum k = L/VS, where A orb = 25+1 
is the number of orbitals in the LLL. We include the 
data for a number of system sizes and rescale the mag- 
netic length Iq by a factor U S/N or -^ to minimize the 
finite size effects. For the model V\ Hamiltonian and 
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FIG. 1. (Color online). The variational energy of the model 
wavefunctions defined by Eqn (2) and (3), against Vi (left 
axis, arbitrary units) and Coulomb Hamiltonian (right axis, 
in units of e 2 /e£s), plotted as a function of momentum. The 
data is generated from system sizes ranging from 6 to 12 elec- 
trons (the inset shows the same plot for the bosonic Laughlin 
state). 



Coulomb Hamiltonian, the dispersion obtained using the 
model wavefunction is in excellent agreement with the 
results from exact diagonalization, both in small k and 
large k regime. Our model wavefunctions compare fa- 
vorably with the exact diagonalization eigenstates, with 
99% overlap for 10 electrons. 

Using the same approach we can construct the 
collective-mode wavefunctions for the entire Rcad-Rezayi 
series. As as an example, we consider an interesting case 
of the Moore- Read state, where in addition to the bosonic 
mode, we also obtain a mode that corresponds to the un- 
paired electron - the neutral fermion (NF) mode [l4[ • On 
the sphere, the root configurations of the two modes are 
given by 

11100100110011 ■ ■ L = 2, 111000110011 ■■■L = 3/2 
11100101010011 ■ ■ L = 3, 111001010011 ■■■£,= 5/2 
11100101010101 ■ ■ -L = 4, 111001010101 ■ ■ -L = 7/2 

; (4) 

The Moore-Read ground-state root partition is given by 
2 electrons in 4 consecutive orbitals [Tpj . Similarly to 
the Laughlin state, any deviation from the uniform back- 
ground density yields the position of the quasihole/quasi- 
electron (not labeled in Eq.©). In Eq.(|4]), the left col- 
umn with an integer angular momentum is the magneto- 
roton mode. The NF mode with a half-integer angu- 
lar momentum is given in the right column. Unique 
model wavefunctions can be constructed by imposing the 
constraint Eq.Q, and in addition a modified constraint 



FIG. 2. (Color online). The variational energy of the model 
wavefunctions for the magneto-roton (MR) mode and the neu- 
tral fermion (NF) mode, evaluated against the 3-body Hamil- 
tonian. The data is generated from system sizes ranging from 
5 to 17 electrons, where the odd number of electrons con- 
tribute to the NF mode, and the even number of electrons 
contribute to the magneto-roton mode. (The inset shows the 
same plot for the bosonic Moore-Read state) 



Eq.© that reads H 3 bCiC2C3\ip^ R ) = 0, where H 3 b is the 
Moore-Read three-body Hamiltonian. Their variational 
energies are plotted in Fig. [2] 

We would like to emphasize for both the Laughlin and 
Moore-Read case, the collective modes enter the multi- 
roton continuum in the long wavelength limit. The con- 
tinuum starts at energy that is double the energy gap of 
the roton minimum. While this makes exact diagonal- 
ization ambiguous, the root partitions give clear physical 
interpretation lor the modes for the entire momentum 
range. Interestingly, the L = 1 state (and L = 1/2 state 
for the neutral fermion mode) vanishes with the set of 
constraints we impose. Thus in the long wavelength limit, 
the collective mode is given by a quadrupole excitation. 
In light of the geometrical picture |l8[ of the FQH effect, 
we identify the magneto-roton mode as a spin-2 "gravi- 
ton", and the NF mode as a spin-| "gravitino", or the 



"supersymmetric partner" of the magneto-roton mode. 

Recently, model wavefunctions for the excitations at 
u = 1/2 filling for bosons were also obtained in a mul- 
ticomponent composite fermion picture [ijl]. We have 
confirmed that the wavefunctions of Rcf. 9( are numer- 
ically identical to ours in finite size systems. The ap- 
proach of Ref . Q , though completely different at the out- 
set, arrives at the wavefunctions that satisfy the same 
clustering properties as ours, which fixes them to be 
unique. The uniqueness property motivates us to inves- 
tigate the SMA wavefunctions obtained from the ground 
state j'i/'o) by the guiding center density modulation, 
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|V>k) = PklV'o), where p k = J2i elkR " is tnc guiding cen- 
ter density operator. The SMA yields excitation energies 
manifestly as a property of the ground state. Though it 
successfully predicts the magneto-roton minimum of the 
collective mode of FQH states at v = 1/m, there is some 
ambiguity in the limit of k — > when the SMA vari- 
ational energy enters the roton-pair continuum, eluding 
comparisons with exact diagonalization. It is thus useful 
to compare the SMA prediction with our model wave- 
functions in the small k limit, since the latter are valid 
in that regime and have transparent physical properties 
given by the root partitions. 

The SMA construction can be adapted to the sphere 
as follows. The ground state on a sphere has the total 
angular momentum L = 0, and the SMA wavefunction 
with total angular momentum L is obtained by boosting 
one electron with orbital angular momentum L. The 
projection into the LLL is equivalent to the projection 
of the boosted single-particle state into the sub-Hilbert 
space of the total spin S. Formally we have 

IV-ir>=E^f+M,M, m J^), (5) 

i 

where i is the electron index, and C r ^f J ' Mm is defined 
by its action on the single electron state C^f^^Jm) = 
C^iLK), where C&f£ m = (m'\Y LM \m) are the 
Clcbsh- Gordon coefficients, and Y LM are the spherical 
harmonics. This is a result of the Wigner-Eckart Theo- 
rem, and due to rotational invariance we can set M = L 
in Eq.©. The dispersion of the SMA wavefunctions is 
plotted in Fig. [3] along with that of our model wavefunc- 
tions. 

For small momenta the variational energies of the two 
classes of wavefunctions agree very well, while the SMA 
mode evidently becomes invalid for momenta larger than 
the magneto-roton minimum. Note that at L = 2,3 the 
SMA wavefunctions only involve the elements of the ba- 
sis squeezed from the same root partition that defines 
our model wavefunctions. Taking the Laughlin 1/3 as an 
example, we now prove the SMA wavefunctions are ac- 
tually identical to ours at L = 2,3. By the product rule 
of the Jack polynomial, we can write 

| Vo) ~ ® J% 2 + Jx 3 ® Jx 4 + \$o) , (6) 

where we suppressed the relative coefficients because 
they are unimportant for the proof. The partitions 
Ai = {10010}, A 2 = {01001001 • • ■ }, A 3 = {10001}, A 4 = 
{10001001 ■ • • }, and |-0o) involves the rest of the squeezed 
basis. It is easy to check that c\Ci C^ L : S L L m .|V>o) = 
0. We thus have 

V lCl c 2 \^ft A ) ~ Vi{0000} ® ( J£ + tJ?J = (7) 

Again, the coefficients are suppressed in Eq.([7]), and t = 
for L = 2. Thus the SMA wavefunctions satisfy ex- 
actly the same constraints as the our model wavefunc- 
tions, which makes them identical. Note that for L > 3 
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FIG. 3. (Color online). The variational energies for the SMA 
model wavefunctions compared to our model wavefunctions 
for Laughlin state at v = 1/3 filling. (The inset shows the 
same comparison for the magneto-roton mode of the Moore- 
Read state) 



the SMA wavefunctions contain unsqueezed basis com- 
ponents with respect to the root partitions used in our 
model wavefunction, and the proof breaks down. 

In conclusion, we demonstrated a numerically efficient 
method of constructing accurate model wavefunctions for 
the collective modes in FQH systems. The wavefunctions 
are identified with the SMA wavefunctions in the long 
wavelength limit. This result reveals a crucial link be- 
tween the "graviton" mode and the SMA mode at long 
wavelengths, which plays an important role in the ge- 
ometrical theory of the FQH effect [la ]. For realistic 
Coulomb or pseudopotential interactions, the "graviton" 
mode decays into multi-roton pairs and appears exper- 
imentally inaccessible. It would be interesting to see if 
the interaction can be tuned in such a way as to expose 
the "graviton" mode at k — > below the roton-pair con- 
tinuum. 

Acknowledgements. We thank B. A. Bernevig for use- 
ful discussions. This work was supported by DOE grant 
DE-SC0002140. 



[1] D. C. Tsui, H. L. Stormer, and A. C. Gossard, Phys. Rev. 

Lett. 48, 1559 (1982). 
[2] R. Laughlin, Phys. Rev. Lett. 50, 1395 (1983). 
[3] The Quantum Hall Effect, 2nd ed., edited by R. E. 

Prange and S. M. Girvin, Springer- Verlag, New York, 

1990. 

[4] G. Moore and N. Read, Nucl. Phys. B 360, 362 (1991). 
[5] N. Read and E. Rezayi, Phys. Rev. B 59, 8084 (1999). 
[6] J. K. Jain, Phys. Rev. Lett. 63, 199 (1989). 



5 



[7] S. M. Girvin, A. H. MacDonald and P. M. Platzman, 
Phys. Rev. Lett. 54, 581 (1985); Phys. Rev. B 33, 2481 
(1986). 

[8] G. J. Sreejith, C. Toke, A. Wojs and J. K. Jain Phys. 

Rev. Lett. 107, 086806 (2011) 
[9] I. D. Rodriguez, A. Serdyniak, M. Hermanns, J. K. 
Slingerland, N. Regnault, arXiv: 1110.0454 
[10] B. A. Bernevig and F. D. M Haldane, Phys. Rev. Lett. 

100, 246802 (2008). 
[11] B. A. Bernevig and F. D. M. Haldane, Phys. Rev. B 77, 
184502 (2008). 

[12] B. A. Bernevig and F. D. M. Haldane, Phys. Rev. Lett. 

102, 066802 (2009). 
[13] N. Read and D. Green, Phys. Rev. B 61, 10267 (2000). 



[14] G. Moller, A. Wojs, and N. R. Cooper, Phys. Rev. Lett. 

107, 036803 (2011); P. Bonderson, A. E. Feiguin, and C. 

Nayak, Phys. Rev. Lett. 106, 186802 (2011). 
[15] F. D. M. Haldane, Phys. Rev. Lett. 51, 605 (1983). 
[16] B. A. Bernevig and N. Regnault, Phys. Rev. Lett. 103, 

206801 (2009). 

[17] The bosonic counterpart has no more than a single par- 
ticle in two consecutive orbitals, yielding a Jack given 
by J^ 2 . which is a function that vanishes as a power of 2 
when two particles approach each other, i.e. the Laughlin 
wavefunction. 

[18] F. D. M. Haldane, Phys. Rev. Lett. 107, 116801 (2011). 
[19] R. Thomale, B. Estienne, N. Regnault and B. A. 
Bernevig, Phys. Rev. B 84, 045127 (2011). 



